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ABSTRACT 

so" 

O According to current models, gamma-ray bursts (GRBs) are produced when 

^ ; the energy carried by a relativistic outflow is dissipated and converted into radi- 

ation. The efficiency of this process, e 7 , is one of the critical factors in any GRB 
■ model. The X-ray afterglow light curves of Swift GRBs show an early stage of 

■^j- ■ flattish decay. This has been interpreted as reflecting energy injection. When 

combined with previous estimates, which have concluded that the kinetic energy 
of the late (> 10 hr) afterglow is comparable to the energy emitted in 7-rays, this 

^ ■ interpretation implies very high values of e 7 , corresponding to > 90% of the ini- 

tial energy being converted into 7-rays. Such a high efficiency is hard to reconcile 
with most models, including in particular the popular internal-shocks model. We 
re-analyze the derivation of the kinetic energy from the afterglow X-ray flux and 
re-examine the resulting estimates of the efficiency. We confirm that, if the flat- 
Q^i 1 tish decay arises from energy injection and the pre-Swift broad-band estimates of 

the kinetic energy are correct, then e 7 > 0.9. We discuss various issues related to 
this result, including an alternative interpretation of the light curve in terms of a 
two-component outflow model, which we apply to the X-ray observations of GRB 
050315. We point out, however, that another interpretation of the flattish decay 

^ ■ — a variable X-ray afterglow efficiency (e.g., due to a time dependence of after- 

glow shock microphysical parameters) — is possible. We also show that direct 
estimates of the kinetic energy from the late X-ray afterglow flux are sensitive 
to the assumed values of the shock microphysical parameters and suggest that 
broad-band afterglow fits might have underestimated the kinetic energy (e.g., by 
overestimating the fraction of electrons that are accelerated to relativistic ener- 
gies). Either one of these possibilities implies a lower 7-ray efficiency, and their 
joint effect could conceivably reduce the estimate of the typical e 7 to a value in 
the range ~ 0.1 — 0.5. 
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Subject headings: gamma-rays: bursts — ISM: jets and outflows — radiation 
mechanisms: nonthermal — X-rays: individual (GRB 050315) 

1. Introduction 

Recent observations by the Swift X-ray telescope have provided new information on the 
early behavior of the X-ray light curve of long-duration (> 2 s) gamma-ray burst (GRB) 
sources. Specifically, it was found (Nousek et al. 2005) that the light curves of these sources 
have a generic shape consisting of three distinct power-law segments oc t~ a : an initial (at 
t < tbrcak.i) with 300 s < tbrcak.i ^ 500 s) very steep decline with time t (with a power- law 
index a,\ in the range 3 < a\ < 5; see also Tagliaferri et al. 2005 and Bartherlmy et al. 
2005); a subsequent (at tbrcak,i < t < tbrcak,2, with 10 3 s < £break,2 ;$ 10 4 s) very shallow decay 
(0.2 < a 2 ^ 0.8); and a final steepening (at t > tbrcak,2) to the canonical power-law behavior 
(1 < ck 3 < 1.5) that was known from pie- Swift observations. 

Nousek et al. (2005) already recognized that these results have direct consequences to 
the question of the 7-ray emission efficiency in GRB sources. This question is important to 
our understanding of the basic prompt-emission mechanism. In the currently accepted inter- 
pretation (e.g., Piran 1999, 2004), the 7-rays originate in a relativistic jet that is launched 
from the vicinity of a newly born neutron star or stellar-mass black hole. In the simplest 
picture, a fraction e 7 of the energy injected at the source is given to the prompt radiation, 
with the remaining fraction (1 — e 7 ) ending up as kinetic energy of ambient gas that is swept 
up by a forward shock and then mostly radiated as early afterglow emission. One attrac- 
tive mechanism for explaining the prompt emission characteristics invokes internal shocks 
that are generated when the ejecta have a nonuniform distribution of Lorentz factors, which 
results in outflowing "shells" colliding with each other at large distances from the source. 

Pre- Swift observations, based on measurements of the 7-ray fluence and of the "late" 
(> 10 hr) afterglow emission, have implied (when interpreted in the context of the basic 
jet model) comparable (and narrowly clustered) values for the radiated 7-ray energy and 
the kinetic energy feeding the afterglow emission, i.e., e 7 ~ 0.5 (e.g., Frail et al. 2001; 
Panaitescu & Kumar 2001a,b, 2002; Berger et al. 2003; Bloom et al. 2003; Yost et al. 2003). 
This result is seemingly problematic for the internal-shocks model, for which an order-of- 
magnitude lower value for e 7 is a more natural expectation (Kobayashi et al. 1997; Daigne 
& Mochkovitch 1998; Kumar 1999; Guetta et al. 2001). This apparent difficulty could in 
principle be overcome if the ejected shells have a highly nonuniform distribution of Lorentz 
factors (e.g., Beloborodov 2000; Kobayashi & Sari 2001). However, to fit the data the shells 
must also satisfy a number of other restrictive conditions, which reduces the attractiveness 
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of this interpretation (see § 5). An alternative proposal was made by Peng et al. (2005), 
who argued that if the jet consists of an ultra-relativistic narrow component (from which the 
prompt emission originates) and a moderately relativistic wide component, with the latter 
having a higher kinetic energy and the former a higher kinetic energy per unit solid angle, 
then the wide component would dominate the late afterglow emission and the 7-ray radiative 
efficiency of the narrow component could be significantly lower than the value of e 7 inferred 
under the assumption of a single-component jet. As explained in Peng et al. (2005), this 
proposal was motivated by observational indications of the presence of two components in 
the late afterglow light curves of several GRB sources and by the predictions of certain GRB 
source models. 

The earliest (steepest) segment of the afterglow light curve is most naturally explained 
as radiation at large angles to our line of sight corresponding either to the prompt emission 
(Kumar & Panaitescu 2000) or to emission from the reverse shock that is driven into the 
ejecta (Kobayashi et al. 2005). This implies that the early afterglow emission is much weaker 
than what would be expected on the basis of an extrapolation from the late-afterglow data. 
This behavior was interpreted by Nousek et al. (2005) as an indication of an even higher 7-ray 
emission efficiency, typically e 7 ~ 0.9. Such a high efficiency could render the internal-shocks 
model untenable. 

Our primary goal is to evaluate the prompt-emission efficiency as accurately and sys- 
tematically as possible on the basis of current data. For this purpose we re-derive in § 2.1 
expressions that explicitly relate e 7 to observable quantities. In particular, we express e 7 in 
terms of the product nf of two parameters, one (k) encapsulating information that could be 
obtained by pie-Swift observations, and the other (/) representing early-time data obtained 
in Swift measurements. In § 2.2 we re-examine the estimates of the kinetic energy during 
the afterglow phase as inferred from the X-ray flux and present a new general formulation, 
correcting errors that have propagated in the literature and have generally led to an un- 
derestimate of the kinetic energy. We then analyze both pie-Swift (§ 2.3) and Swift (§ 2.4) 
data in a uniform manner in the context of this formalism. We argue that the kinetic energy 
estimates remain subject to considerable uncertainties. In particular, while the simple anal- 
ysis of a large number of bursts (using "typical" values for the microphysical parameters) 
yields rather high values for the kinetic energy and hence a low inferred 7-ray efficiency, a 
mult i wavelength analysis (which determines the microphysical parameters from the fit to the 
data) of a small subset suggests that the kinetic energy is lower and hence the inferred value 
of e 7 is higher. Some caveats to this analysis are considered in § 2.5. 

The suggestion that the Swift observations imply a higher value of e 7 in comparison with 
the pie-Swift results is based on the interpretation of the flatfish segment of the X-ray light 
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curve as reflecting an increase in the kinetic energy of the forward shock during the early 
stages of the afterglow. In this picture, the kinetic energy just after the prompt-emission 
phase was significantly lower than the kinetic energy estimated from the later stages of the 
afterglow (the pie-Swift results). One conceivable way of avoiding the need to increase the 
estimate of e 7 in light of the new Swift observations is through a time evolution (specifically, 
an increase with t) of the X-ray afterglow emission efficiency ex- Such a behavior could in 
principle also account for the flatfish segment of the light curve and eliminate the need to 
invoke an increase in the shock kinetic energy. If, in addition, the value of the afterglow 
kinetic energy at late times were underestimated by the broad-band fits to pie- Swift GRB 
afterglows, which could be the case if only a fraction £ e < 1 of the electrons were accelerated 
to relativistic energies in the afterglow shock (see Eichler & Waxman 2005), then the typical 
afterglow efficiency would be further reduced (to a value as low as e 7 ~ 0.1 if £ e ~ 0.1), 
which might reconcile the new Swift data with the comparatively low efficiencies expected 
in the internal-shocks model. We discuss these issues in § 3. 

A two-component jet model with the characteristics required for reducing the inferred 
7-ray emission efficiency is evidently disfavored by the Swift data, but a model of this type 
with different parameters could provide an alternative interpretation of the flatfish shape of 
the light curve between tbreak.i and £break,2- We elaborate on these matters in § 4, where we 
also present a tentative fit to the X-ray light curve of the Swift source GRB 050315 in the 
context of this scenario. 

Our conclusions on the physical implications of the early X-ray light-curve observations 
of GRB sources are presented in § 5. 

2. Estimating the Gamma- Ray Efficiency 

The observed isotropic-equivalent luminosity can generally be expressed as L iso = eE iso /T, 
where E iso is the isotropic-equivalent energy of the relevant component, T is the duration of 
the relevant emission, and e is the overall efficiency. This efficiency is a product of several 
factors: e = edi S e e eradeobs, where a fraction edis of the total energy is dissipated into internal 
energy, a fraction e e of the internal energy goes into electrons (or positrons) and can in 
principle be radiated away, a fraction e ra d of the electron energy is actually radiated, and 
a fraction e Q b s of the radiated energy falls within the observed range of photon energies. A 
fraction e 7 = edis,GRBe e ,GRBerad,GRB of the total original (isotropic equivalent) energy -Ei SO ,o is 
radiated away during the prompt emission, -E 7i i SO = e 7 -Eiso,o, while the remaining (kinetic) 
energy in the 7-ray emitting component of the outflow, -Ek,iso,o = (1 — e 7 )i?i SOi o, is transferred 
to the afterglow shock at tdcc- In the prompt GRB emission Cgrb = ^obs^RB- In the 
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afterglow edi S ,x ~ 1, and therefore ex ~ e ei xe ra d,xeobs,x, where we concentrate on the X-ray 
afterglow. 

Using the above expression for e 7 , it becomes clearer why it is difficult for it to assume 
very high (> 0.9) values. Whereas e radiG RB ~ 1 is possible, and even expected if the elec- 
trons cool significantly over a dynamical time (as is typically expected in the internal-shocks 
model), e 7 > 0.9 requires in addition that both edi S ,GRB > 0.9 and e ej GRB > 0.9. It is difficult 
to achieve edi S ,GRB > 0.9 (i.e., dissipate more than 90% of the total energy) in most models 
for the prompt emission, and in particular in the internal-shocks model. Furthermore, it is 
not trivial to produce e ej GRB > 0.9 (i.e. more than 90% of the postshock energy going to 
electrons) in the internal-shocks model, where the electrons are believed to be accelerated 
in a shock propagating into a magnetized proton-electron plasma. In particular, this would 
require that less than 10% of the postshock energy goes into the protons and the magnetic 
field, which is difficult since the protons carry most of the energy (in kinetic form) and the 
magnetic field likely carries a comparable energy in the upstream fluid ahead of the shock. 
It might conceivably be possible if less than 10% of the energy is in the magnetic field and 
the protons can somehow transfer their energy to the electrons, which radiate it away. 



2.1. Relationship to Observed Quantities 

The X-ray afterglow isotropic-equivalent luminosity, Lx,i SO , can be expressed in terms of 
the X-ray afterglow flux, Fx, 

L x ,iso(t) = 4tt4(1 + zf'^Fxit) , (1) 

if F u oc v~Pt~ a in the relevant ranges in frequency and time. 1 The efficiency of the X-ray 
afterglow emission is defined as 

ex(i) - (2) 
where -Ek,iso is the isotropic-equivalent kinetic energy in the afterglow shock. 
Using the relation E isofi = -E 7iiso /e 7 = E Kisofi /(l - e 7 ), we obtain 

rpobs Z?obs p /, \ 

fc 7 fc obs,GRB 7,iso ^ 7,iso j J - / k,iso^' / */ 

(1 — e 7 ) -E^isc^O -^k,iso(£*) -E^iso.O 



1 For the more general case, see eq. [1] of Nousek et al. (2005) and the discussion thereafter. Here 
Lx.iso(t) = J^ 2 dv L„.i so (t), where both v and t are measured in the cosmological frame of the GRB, whereas 
Fx{t) — J^ 2 dvF v (t), where both v and t are measured in the observer frame. 
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where i* can be chosen as a time at which it is convenient to estimate -Ek,iso, and we shall 
use i* = 10 hr (since it is widely used in the literature and is typically > tbrcak,2, the end of 
the flattish segment of the X-ray light curve). 

An important question is how to estimate / and re from observations. The most straight- 
forward way of estimating / = Fk ; i so (10 hr)/Fk ; i so ,o, which has been used by Nousek et al. 
(2005), is by evaluating the value of the X-ray flux decrement at tdcc relative to the extrapo- 
lation to tdec of the late-time (t > tbreak,2) flux, and translating this flux ratio into an energy 
ratio according to the standard afterglow theory. When estimating / in this way we assume 
that 

e x ~(l + F)- 1 e e K/z/x) {p - 2)/2 , (4) 

where v m and v c are the characteristic synchrotron frequency and cooling break frequency, 
respectively (Sari et al. 1998; Granot & Sari 2002) and Y is the Compton y-parameter. 
This result can also be obtained under the assumptions of standard afterglow theory, as 
follows. The overall X-ray afterglow efficiency can be written as ex ~ e e e ra de bs, where 
6 rad » min[l> m M)^ 2 )/ 2 ] and e obs » (1 + F) _1 max[(z/ m /z/ x ) (p ~ 2)/2 , (*V^x) {p ~ 2)/2 ]; the 
factor of (1 + Y)~ l is the fraction of the radiated energy in the synchrotron component, 
and it is present because the synchrotron self-Compton (SSC) component typically does not 
contribute significantly to -Fx (10 hr) but may still dominate the total radiated luminosity. 
The factor (1 + Y)^ 1 is generally assumed to be small, consistent with the usual inference 
that the magnetic-to-thermal energy ratio in the emission region, e#, is smaller than e e (see 
eq. [12]). 2 However, one should keep in mind that a different time dependence of ex (which 
might occur under less standard assumptions) would modify the value of / accordingly. The 
value of e bs,GRB can be estimated by extrapolating the observed part of the spectrum and 
modeling additional spectral components that might carry considerable energy. 

To estimate re = F°j! o /F kiiso (10 hr) we calculate E° h * = j^-nd\ (1 + z)~ l directly from 
the observed 7-ray fluence, / 7 , and the measured redshift, z. The denominator, F ki i SO (10 hr), 
has been estimated following Freedman & Waxman (2001) and Kumar (2000) (see also 
Lloyd- Ronning & Zhang 2004) from Fx (10 hr) using standard afterglow theory. We recon- 
sider this calculation in the next section. 



2 At early times (less than about a day) inverse-Compton radiation is important in cooling the electrons 
that are emitting synchrotron X-rays. Later the X-ray emitting electrons are within the Klcin-Nishina cutoff 
and are not cooled by inverse-Compton radiation (see Fan & Piran 2006). 



-7- 



2.2. Estimates of the kinetic energy from the X-ray afterglow observations 

Using the results of Granot & Sari (2002) we find that, if vx > max(i/ m ,j/ c ), we can 
estimate the kinetic energy from the observed X-ray flux. We use 

P W-QOy in52 g(g) /-, , V \4/(p+2) 4(l-p)/(p+2) (2-p)/(p+2) r 4/(p+2) .(3p-2)/(p+2) /rx 

^k,isoW — y-^ x iu ^ 2) 1-1 B ~ 2 x ' 46 10hr 

where Lx.iso = -^x,46l 46 erg s -1 is the isotropic-equivalent X-ray luminosity in the range 
2 — 10 keV at a time 10ti h r hr, both measured in the cosmological frame of the GRB 
[corresponding to an observed time of t = 10(1 + z)t±ohr hr and spectral range 2/(1 + z) — 
10/(1 + z) keV], e e ,_i = e c /0.1, e B ,_ 2 = e B /0.01, and 



g(p) 



(p- l) J, - 1 exp(5.89p- 12.7) 
( 5 (p-2)/2 - i)(p - 2)P" 2 (p - 0.98) 



4/(p+2) 



(6) 



In equation (5), Lx,iso is evaluated through Fx using equation (1), and it is assumed that 
Fx is dominated by the synchrotron component. If there is a significant contribution to 
Fx(t*) from the SSC component then equation (5) would overestimate -Ek,iso(£*)> but the 
correct estimate could still be obtained if only the synchrotron contribution to Fx(t*) is used 
(although in practice it might be difficult and somewhat model-dependent to disentangle 
the synchrotron and SSC components). For our fiducial parameters (1 + Y) m 3.7, so the 
numerical coefficient in equation (5) is 3.2 x 10 53 erg for p = 2.2 and 5.8 x 10 53 erg for 
p = 2.5. This is a factor ~ 30 — 60 higher than the numerical coefficient in equation (7) of 
Lloyd- Ronning & Zhang (2004). Part of this difference (a factor of ~ 3 — 3.5) reflects the 
fact that these authors did not take into account the SSC contribution, which reduces the 
flux in the relevant power-law segment of the spectrum by a factor of (1 + Y) (Sari & Esin 
2001; Granot & Sari 2002), 3 but this does not account for most of the discrepancy. Most of 
the difference is basically a result of a higher (by a factor of ~ 38 for p = 2.2) value of v m 
that Lloyd- Ronning & Zhang (2004) use in their equation (2), given that their expressions 
for v c and F umax (their eqs. [3] and [4]) are very similar to those in Granot & Sari (2002) 
and that F„ >max(Wc) = F^yj 2 ^^ 2 ^/ 2 oc u { ^ 1)/2 . Our numerical coefficient is lower 
than that in equation (4) of Freedman & Waxman (2001) by a factor of ~ 3 and ~ 8 (or 
~ 12 and ~ 31 if SSC is not taken into account) for p = 2.2 and p = 2.5, respectively. The 
difference in the numerical coefficient between Lloyd- Ronning & Zhang (2004) and Freedman 
& Waxman (2001) is by a factor of ~ 100 and ~ 475 for p = 2.2 and p = 2.5, respectively. 



3 This is assuming that the SSC component does not contribute considerably to the observed X-ray flux, 
which is typically the case at t* = 10 hr. 
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Using equations (2) and (5), we can express the efficiency as 



e x (t) = 3.5 x 10" 



i , V^- 1 ^ -1 t (f- 2 )/ 4 P (p-2)/4 ,-3(p-2)/4 
1 ~t~ 1 ) t e,-l t B-2 "^k.iso^ hohr 



(7) 



- QQyin^ ^( 2 - 2 ) (i , V \-4/(p+2) 4(p-l)/(p+2) (p-2)/(p+2) f (p-2)/(p+2) ,-2(p-2)/(p+2) 
_ ~ l ~ ' e - 1 S -2 ^X,46 C 10hr ' 

where -Ek,iso,52 = -^k,iso(^)/(10 52 erg). For our fiducial values (and for Lx,46 rather than 
£k,iso,52), the numerical coefficient in equation (7) is 1.1 x 1CT 3 for p = 2.2 and 6.3 x 1CT 4 
for p = 2.5. 

Equations (5)-(7) are valid for p > 2, but they can be easily generalized to p < 2 by sub- 
stituting e e -> e e (p-l)/(p-2), where e e = e e 7 min /(7 e ), (7 e ) = / d'j e (dN/d'j e )'j e [J d'j e (dN/d'j e )] 
is the average electron Lorentz factor, and the electron energy distribution is a power law of 
index p at low energies above 7 m i n . If the electron energy distribution is a single power law, 
dN/dj e oc % p for 7 min < 7 e < 7 max , then 



Q (7e) fP~ l \ 1 - (7max/7min) 2 P 



7min \P 2/ 1 (7max/7min) ^ 



»(p-l)/(p-2) 


p>2, 


ln(7max/ 7min) 


p = 2, 


^(7max/7min) 2 - p (p-l)/(2-p) 


1 < p < 2 , 


(7max/ 7min)/ ln(7 max / 7min) 


p = 2, 


~ (7max/7min)(l ~ p) / (2 ~ p) 


p < 1 . 




(8) 



This results in slightly modified equations: 

^,iso(0 = 1-19 X 1Q 52^ (1 +r) 4/( P+ 2) e -4(l- P )/( P+ 2) e g- P )/( P+ 2 ^ (g) 



3.01 x 10" 
3.03 x 10" 



52 


0(p) n 


0(2.2) v 




0(p) = 


-2 


" 0(p) " 




[0(2.2). 


-3 0(2-2) 




0(P) ' 



(10) 



(p - 2) exp(5.89p- 12.7) 1 4/(p+2) 
(5(p-2)/2 _ i)(p_o.98) J 

" P+2)/4 (i + ^)- 1 C-\4r 2 )/4 «^or 2)/4 (id 

| , v^- 4 /(P+2)= 4 (P- 1 )/(P+ 2 )>-2)/(P+ 2 ) I-(P-2)/(P+2),-2(p-2)/(p+2) 
/ fc e,-l fc _B,-2 -^X,46 hohr 



Note that the numerical coefficient was calculated in Granot & Sari (2002) only for p = 
2.2, 2.5, 3 and interpolated between these values. Extrapolating that formula to p < 2 could 
potentially be very inaccurate. 
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For simplicity we use the expression for Y that is valid in the fast-cooling regime, 

Y £t/£i,>>1 , (12) 

I e e /e B e e /e B < 1 

which is still reasonable at 10 hr, even if it is slightly after the transition to slow cooling. 
More generally, e e /es should be multiplied by e ra d ~ min [l, (f m /V c )( p ~ 2 )/ 2 ] , where for p < 2 
and u c < z/ max we have e ra( j ~ 1. 

Before applying these relations to observed bursts we remark on a common misconcep- 
tion concerning the dependence of -Ek.iso that is inferred from Fx on and e e . It has been 
argued that £k,iso is very insensitive to the exact value of €b (e.g., Freedman & Waxman 2001; 
Piran et al. 2001). This follows from the observation (see eq. [5]) that for vx > max(z/ m , u c ), 

E Klso OC (l + y)4/(p+2) e j4(p-l)/(p+2) e -(p-2)/(p+2)^ suggestg thftt ^ depends yery weakly 

on 6b- However, this holds only in the limit where (1 + Y) ~ 1, which corresponds to Y <C 1 
and €b ^> e e , whereas observations suggest that the opposite limit typically applies, e# <C e e , 
in which case (1 +Y) F fa (eg/ee) 1 / 2 and i?k,iso oc e B p ^ p+2 \ This is a significantly stronger 
dependence on e^. (Note that the inferred value of varies by ~ 2 orders of magnitude 
among different afterglows, from ~ 10~ 3 to ~ 0.1, corresponding to a change of more than 
an order of magnitude in the estimated value of £k,iso-) The dependence of -Ek,iso on e e is 
stronger, oc e ~ 4 ^'~ 1 ^^ p+2 ^ (i. e . ; slightly steeper than an inverse linear relation) in the limit 
£b ^ e e , but only oc e e 2ly2p 3 ^ p+2 ) in the more relevant limit of eb <C e e . In the latter case 
the dependence of -Ek,iso on e B is rather similar to its dependence on e e . 4 However, e e appears 
to vary much less than e# among different afterglows, only covering a range of about one 
order of magnitude (between ~ 10~ L5 and ~ 10 -0 ' 5 ), which corresponds to a variation in 
-E^iso by a factor of ~ 4. It is also worth noting that the expression for -E^iso has some 
(nontrivial) dependence on the value of p (see eq. [6]). 



2.3. Pre-Swift GRBs 

Table 1 shows the estimated values of k deduced from the observational properties of 
17 pie-Swift GRBs with known redshifts, using the samples of Berger et al. (2003) and 
Bloom et al. (2003). We provide the values .E'k.iso.iohr = -Ek,iso(10 hr) and k for our fiducial 
parameter values (e e = 0.1, es = 0.01, and p = 2.2). The value of -E^iso.iohr and therefore 



4 Quite often Y ~ 1 — 2 is inferred, in which case neither of the asymptotic limits F <C 1 and Y 3> 1 is 
applicable and the dependence of -Ekaso on and e e does not have a power-law form but rather the more 
complex form given by eqs. (5) and (12). 
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of k depends on the values of the microphysical parameters (e e , e#, and p) that are not 
well known. Therefore we also calculate -Ek,iso,iohr and k (for those GRBs for which this is 
possible) using the values of the microphysical parameters inferred from the fits to the broad- 
band afterglow data that were performed by Panaitescu & Kumar (2002, denoted by PK02, 
using Table 2 therein) and by Yost et al. (2003, denoted by Y03, using Table 1 therein). 
Finally, we compare the values we obtain for -Ek,iso,iohr using equation (5) to those obtained 
for £k,iso(l day) = £k,iso,id by Y03 and those obtained for E kiiso (10 hr) 0.5-E k ,i S o,o by PK02. 
Also shown in Table 1 are the corresponding values of k, including K ld = E 1: - lso /E^ old . 

When using the fiducial parameters e e = 0.1, €b = 0.01, and p = 2.2 our estimates for 
-E'k.iso.iohr are significantly larger than the estimates of Lloyd- Ronning & Zhang (2004), who 
use the same values. This can be traced to the difference in the numerical factor that appears 
in equation (5). These relatively large values of -Ek,iso,iohr lead to a typical value of k around 
0.1 — 0.2, for which the 7-ray efficiency problem would not be very severe. Similar results 
were obtained by Fan & Piran (2006), whose estimates for £k,iso,iohr are within a factor of 2 
of those presented here. 

The situation is different when we use the values of the microphysical parameters from 
the pie-Swift afterglow fits. In these cases the values of -Ek,iso,iohr are typically lower, resulting 
in higher estimates for k. The estimates of the kinetic energy from the PK02 fits to the 
afterglow data (rather than from using eq. [5]) are generally the lowest. 5 The corresponding 
values of k are close to unity. These results reflect the pie-Swift inference that there exists a 
rough equality between the isotropic-equivalent 7-ray and late (10 hr) kinetic energies. The 
comparison between the detailed calculations and those based on equation (5) may suggest 
that we might be better off adopting different fiducial parameters (e.g., e e = 0.3, e# = 0.08, 
p = 2.2). We also note that the values of -Ek,iso,iohr obtained by substituting the values of 
the microphysical parameters from the afterglow fits into equation (5) are generally higher 
(corresponding to lower values of k) compared to those derived directly by those fits. The 
"typical" values of the microphysical parameters inferred from the fits are roughly e e ~ 0.3, 
e B ~ 0.03, p « 2.2. 

The values of k obtained here are crucial to the overall estimate of e 7 and to the origin 
of the "high efficiency problem." Considering Table 1, one should proceed with care in view 
of the large dispersion in the different estimates of -Ek,iso,iohr- An example is a factor of 4.7 
between the independent estimates of PK02 and of Y03 for GRB 000926 (and a factor of 
2.5 in the opposite direction for GRB 970508). The dispersion is even greater between the 



5 In PK02 the values of E^- 1SO (10 hr) are typically a factor of ~ 2 smaller than i? k ,iso,o due to radiative 
losses at early times (A. Panaitescu, personal communication; no energy injection is assumed in that work). 



estimates of -Ek,iso from the afterglow fits and the values obtained using equation (5) for the 
same values of the microphysical parameters — a factor of 18 for Y03 (and 11 for PK02) for 
GRB 000926 . 



2.4. 



Swift GRBs 



The new result found by Swift is the appearance of a rapid decay followed by a shallow 
decline phase in the X-ray afterglow. For seven out of the ten Swift GRBs considered here 
there was a clear observation of the two breaks in the light curve, at tbreak.i and tbreak,2 
(the beginning and end, respectively, of the flattish segment of the X-ray light curve). If 
we interpret the shallow decline as arising from an additional injection of energy into the 
afterglow shock (see § 1) then the ratio of the X-ray fluxes at the end and at the beginning 
of this phase can be used to estimate / (the corresponding ratio of the kinetic energies) for 
these bursts (Nousek et al. 2005). A lower limit, / min , is obtained using the flux decrement 
at ^break.i relative to the extrapolation back to that time of the late-time (t > £break,2) flux. 
An approximate upper limit, /max, is obtained by assuming that the flat part of the emission 
from the forward shock starts at tdcc ~ ^grb and is simply buried underneath the tail of 
the prompt emission at t < t break, i- While formally f min < f < /max, it is reasonable, in 
the context of the basic jet model, that the assumption made to calculate / ma x is basically 
applicable, so that / ~ /max- Under this assumption one infers / > 10, and in some cases 
even a much larger value of / (~ 10 2 — 10 3 ). 

The results for / need to be combined with an estimate of k. Table 2 shows the values 
of -Ek,iso,iohr and k for the ten Swift GRBs with known redshifts from the Nousek et al. 
(2005) sample, estimated using equation (5) with e e = 0.1 and ee = 0.01. The value of p 
was derived using the measured spectral slope in the X-ray band (attempting to fit it into 
the range 2 < p < 3 if allowed within the measurement errors). Lacking any broad-band fits 
to Swift bursts, this is the best direct evidence that we have so far from these data. The 
resulting values of k are similar to those from the pie-Swift era (see Table 1). With the 
exception of GRB 050401, for which k = 0.41, and a few other bursts for which we only have 
upper limits, the inferred values of k are less than 0.1. If this is the correct value of k then, 
using 



(see eq. [3]), we find that with / ~ 10 the overall 7-ray efficiency is not larger than ~ 0.5 
(assuming e bs,GRB ~ 1). A similar conclusion was reached by Fan & Piran (2006). 

One may question this conclusion in view of the fact that, in pre-Swift bursts, broad- 
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band analyses of the afterglow data generally resulted in a significantly lower values of 
-Ek,iso,iohr, and correspondingly higher values of re, compared to those obtained from equation 
(5) with the same fiducial values of the microphysical parameters (e e = 0.1 and e# = 0.01; 
see Table 1). Furthermore, the choice of the fiducial values of the microphysical parameters 
is somewhat arbitrary and it affects the resulting values of -Ek,iso,iohr and re. It is reasonable 
to expect that the values of the microphysical parameters that would have been inferred 
from a broad-band fit to the afterglow data of the Swift bursts would have led to higher 
estimates of -Ek,iso,iohr and re that were similar to those derived for the pie-Swift GRBs. The 
latter values, however, vary among the different estimates, from as high as ~ 5 — 8 to as 
low as ~ 0.1 — 0.3 (see Table 1). In light of this, one may adopt a "typical" value of re ~ 1, 
keeping in mind that there is an uncertainty of almost an order of magnitude in this value. 

The adoption of this higher value of re (~ 1) for the Swift GRBs, similar to the values 
inferred from broad-band modeling of pie-Swift sources, together with the interpretation 
of the shallow decay phase as arising from energy injection (and hence / ~ 10) leads to 
the conclusion that typically e 7 ~ 0.9, and in some cases e 7 is even as high as ~ 99% (or, 
equivalently, 1 — e 7 e bs,GRB/«/ is as low as ~ 10~ 3 — 10~ 2 ). Such a high 7-ray efficiency 
would be extremely hard to produce in the internal shocks model (see § 1). 

If, on the other hand, re ~ 0.1 and there is energy injection (i.e. / ~ 10), or if re ~ 1 
and there is no energy injection (i.e. / = 1; see § 3), then this would imply a significantly 
lower typical 7-ray efficiency, e 7 ~ 0.5, although in some cases e 7 would still be as high as 
~ 90% (or, equivalently, 1 — e 7 ~ e bs,GRB/ K / would still be as low as ~ 10~ 2 — 10 _1 ). Even 
the latter, more moderate, requirements on the 7-ray efficiency are not easily satisfied in the 
internal-shocks model, athough they might possibly still be accommodated in this scenario 
(Kobayashi et al. 1997; Kumar 1999; Guetta et al. 2001). Finally, if re ~ 0.1 and the shallow 
decline does not arise from energy injection (i.e., / = 1) but, say, from varying afterglow 
efficiency (as discussed in § 3), then the 7-ray efficiency would typically be e 7 ~ 0.1. In the 
latter case the results are fully consistent with the predictions of the internal-shocks model. 



2.5. Some Caveats 

The discussion so far was relevant to the power-law segment of the spectrum where 
ux > max(z/ m , u c ). If, instead, u m < z/ x < u c , we have only a lower limit on the value of E^ iso 
from this consideration (i.e., from eq. [5]) that corresponds to an upper limit on the value 
of re, the true value being smaller than this upper bound by a factor [f c (10 hr) / ux] 2 ^ p+2 ^ 
(which, however, is not typically expected to be very large). 
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It is possible that there is a contribution to Fx(tdec) from a component of the outflow that 
is not along the line of sight, or for some other reason did not contribute to the observed 7-ray 
emission. In this case equations (3) and (5) would provide a lower limit on e 7 , rather than 
directly determine its value. If both the prompt 7-ray emission and -Fx(tdec) were dominated 
by emission from angles 9 > 1/T relative to the line of sight (where T is the Lorentz factor 
of the outflow), then again equations (3) and (5) would provide a lower limit on e 7 , since 
the beaming of radiation away from the line of sight is expected to be either comparable or 
somewhat smaller during the afterglow emission at t dcc compared to the prompt GRB. 

It is also possible that some fraction e* of the original energy -Eiso.o ended up in a totally 
different form, i.e., was not radiated during the prompt emission and did not end up in the 
kinetic energy of the afterglow shock. This could occur, for example, if along some directions 
a forward shock is not formed (or at least not formed efficiently) for a very pure Poynting- 
flux outflow. In such a case some of the energy (potentially even most of the energy) could 
be carried out to very large distances (in principle out to infinity) in electromagnetic form 
(such as low-frequency electromagnetic waves). 6 Alternatively, a good fraction of the energy 
might be carried away in high-energy cosmic rays and neutrinos (Waxman 1995; Waxman 
& Bahcall 1997) and thus would not contribute to the kinetic energy of the afterglow shock. 
In this case we have -Ek,iso,o — (1 — e 7 — e*)i?iso,o and we need to make the substitution 
(1 — e 7 ) — > (1 — e 7 — e*) everywhere, so the estimate (13) for e 7 (assuming e* = 0) should be 
multiplied by (1 — e*), 



This means that e 7 < 1 — e* (even for nf ^> e bs,GRB), and therefore e* > 0.5 would imply 
e 7 ~ 0-5 • Thus, even when one infers e 7 / (1 — e 7 ) 1 and hence 1 — e 7 <C 1 under the usual 
assumption that e* = (or at least e» < l-e 7 ), the intrinsic 7-ray efficiency might still be 
significantly smaller, and is e 7 < 0.5 for e* > 0.5. Note that, in order for GRBs to produce 
the highest-energy cosmic rays, their energy should be comparable to that of the prompt 
7-rays (Waxman 1995, 2004), i.e., e 7 < e* (the inequality arising since there might be other 
forms of energy that escape the prompt emission site without being directly detected), and 
therefore e 7 < 0.5. 

One should, however, keep in mind that high-energy cosmic rays and neutrinos must 
tap the same dissipated energy that also powers the prompt 7-ray emission. Therefore, for 
the same observed energy in 7-rays and inferred kinetic energy in the afterglow, in addition 
to increasing the required total initial energy E iso fi by a factor of (1 — e*) -1 , these particles 



6 Both the formation of a forward shock and the ability of energy to escape to infinity in electromagnetic 
form have not yet been fully worked out (e.g., Melatos & Melrose 1996; Smolsky & Usov 2000). 
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would also increase the required dissipated energy edis,GRB-Eiso,o and the fraction edi S ,GRB 
of the dissipated energy that ends up in 7-rays. 7 In other words, edi S ,GRB > 67 + 6* = 

(-^7,iso ■E*,iso)/ (-^k,iso,0 -^7,iso -^*,iso) ^ ^y,iso/ (-^k,iso,0 -^7,iso)) where E^^q £*-Fi SO) o 

is the (isotropic-equivalent) energy in cosmic rays and neutrinos, and i?k,iso,o an d -E 7iiso are 
determined (at least in principle) by observations. Clearly, E^^so increases the lower limit on 
edis.GRB- Nevertheless, the fact that E*^ so may reduce e 7 to < 0.5 even for nf >> 1 makes it 
possible to have e ejG RB ^ 0.5, which should be easier to accommodate for shock acceleration 
in a proton-electron plasma. Still, dissipating and getting rid of almost all of the energy 
(through radiation, cosmic rays, neutrinos, etc.) and leaving only a small fraction of the 
original energy in the kinetic energy of the forward shock, as is required for nf >> 1, is not 
an easy task for any model of the prompt emission. 



3. Efficiency of the X-Ray Afterglow Emission 

As noted in § 1, one of the new features discovered by Swift is the early shallow decline 
phase: F x oc t~ a with 0.2 < a < 0.8 for t b rcak,i < t < t b reak,2- During this phase tF x (t) 
increases with time. Using the definition of exit) (eq. [2]) and the relation between Lx,iso 
and F x (eq. [1]), we find that 

ex{t ^ {t) = 4.4(1 (is) 

is constant in time. Note that if v and t in the expression for Fx(t) were referred to the GRB 
rest frame rather than to the observer frame then the factor (1 + z)P~ a ~ l on the right-hand 
side of equation (15) would be eliminated, and with it any potential (weak) time dependence 
resulting from a possible temporal variation of a or f3. 

Now, if the observed frequencies satisfy vx > max(z/ m , v c ), and p > 2, then equation (4) 
is applicable. Furthermore, if the afterglow shock evolves according to the adiabatic self- 
similar solution of Blandford & McKee (1976) and the fractions of the postshock internal 
energy in electrons (e e ) and in magnetic field (e#) are constant in time, then v m oc t~ 3 / 2 . 
Under these circumstances ex oc t~ 3 ( p ~ 2 )/ 4 decreases (slowly) with time. As can be seen 
from equation (15), €x(t) E^i so (t) oc tFx(t). Therefore, the observed rise in tFx(t) implies a 
similar rise in ex{t) E^i SO (t) . Given the expected decrease of ex{t) with t for p > 2, the slowly 
decaying portion of the light curve has been attributed by several researchers to an increase 



7 In the case of a very pure Poynting flux, the escaping energy E* t i so is the fraction that did not dissipate. 
Therefore, while £?i so ,o increases by a factor of (1 — e») _1 , the dissipated energy £dis,GRB-Eiso,o remains 
unchanged (assuming other efficiencies do not change), and thus edi s ,GRB decreases by a factor of (1 — e*) _1 . 
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in -E'k,iso(0) i- e -' to some sort of energy injection into the forward shock (e.g., Nousek et al. 
2005; Panaitescu et al. 2005; Zhang et al. 2005; Granot & Kumar 2006). 

It is, however, conceivable that the rise in tF x (t) corresponds, at least in part, to an 
increase of ex(t) with time. One way in which this could be brought about is if p were < 2 
(assuming e e and e B remain constant). 8 In this case N(j e ) oc % p for 7 min < 7 e < 7 max and 
ex ~ (1 + Y)^ 1 e e (u miiX / uxY P ~ 2 ^ 2 [the same as eq. (4) for p < 2, but with z/ max replacing z/ m ], 
where z/ max = ^ syn (7 max ) oc 7-B7m ax oc 7 4 Pext oc t~ 3 / 2 [where B is the comoving magnetic 
field amplitude and p cxt is the external density; the same as the scaling of v m = ^ S yn(7mm) 
for p > 2] , so ex oc 

t 3(2- P )/4_ Thig time 

dependence is the same as for p > 2, but in this case 
ex increases with time whereas for p > 2 it decreases. Similarly, F x oc t - ^" 2 )/ 4 , just as for 
p > 2, except that for p < 2 this corresponds to a decay rate flatter than t~ l . One possible 
difference between the two cases is that for p < 2 and v c < z/ max radiative losses are not 
always negligible since most of the energy in the electrons is at 7 e ~ 7 max and is therefore 
radiated away. Thus, unless e e < 1, radiative losses could be significant and would tend to 
steepen the flux decay rate and make it harder to achieve a flatfish light curve. We also 
note that the X-ray spectral slope for p < 2 [assuming z/x > max(z/ c , u m )] is /3x = p/2 < 1, 
which in many cases is inconsistent with the observed value (Nousek et al. 2005), so this 
explanation of the shallow decay of F x might only apply to a subset of the sources (see Fig. 
8 of Nousek et al. 2005). 

An alternative possibility for ex(t) to increase with t is for either one (or both) of the 
microphysical parameters e e and e B to increase with time. Using the dependence of v m 
on these parameters (e.g., Sari et al. 1998), we find that, for p > 2, ex oc ^e^ 2 ^ 4 
when e e < and ex oc e^" 3 ^ 2 ^^ 4 when e e ^> e B . This applies when parameterizing 
in terms of -Ek.iso (which is not measured directly), whereas a parameterization in terms 
of Lx,iso (which is measured directly) yields ex oc e 4 ( p ~ 1 )/( p+2 ) e ^'~ 2 )/( p+2 ) when e e e# and 
e x oc e 2 ( 2 P- 3 )/(p+ 2 ) e P/(p+ 2 ) w h en e e > e B (see eq. [2]). Table 3 of Nousek et al. (2005) provides 
the values of Act — the change in the temporal decay index across the break in the light 
curve at tb re ak,2 — for nine Swift GRBs in which it could be measured reliably, and shows 
that typically 0.5 < Aa < 1. In our context, if e e oc t ae and e# oc t aB at t < tbrcak,2, then 
attributing the flatfish decay phase to a growth of e e and/or e# with time requires (in the limit 
e e > e B ) that Aa = 2(2p - 3)/(p + 2)a e + a B p/(p + 2). For 2 < p < 3, 0.5 < p/(p + 2) < 0.6 
and 0.5 < 2{2p - 3)/(p + 2) < 1.2. Therefore, for p w 2, a e + a B ~ 2Aa ~ 1 - 2 and 



8 It is in principle possible that ex could increase with time on account of its dependence on Y (see eq. 
[4) even if e e and es remained constant and p were > 2, given that Y decreases with time in the slow- 
cooling regime (y m < v c ) that is relevant for p > 2. However, one can show that, to be relevant during the 
early-afterglow phase, this would require unrealistically high values of p. 
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for Aa ~ 0.5, a roughly linear growth with time of either e e or €b (or of their product) is 
required. For p ~ 2.6 and Aa ~ 1, a linear growth of e e and a constant e# would work (i.e. 
a e ~ 1 and = 0). If Y ~ 1 (rather than 7 > 1 or F C 1), the dependence of e x on 
e e and e# is no longer a simple power law, requiring a similarly nontrivial dependence of e e 
and/or e# on the observed time t (insofar as Fx is indeed a pure power law in t during the 
flatfish decay phase). A physical scenario will need to account both for this behavior and 
for why the time evolution of the microphysical parameters effectively terminates at t D reak,2- 

The magnetic-energy parameter e# could reflect either the structure of the ambient mag- 
netic field or postshock field- amplification processes. In the former case an increase of e# 
with time could be caused by an increase of the magnetization parameter a = -Bg Xt /47rp cxt c 2 
of the ambient gas with distance from the source, which might occur in certain GRB progen- 
itor models (e.g., Konigl & Granot 2002). In the latter case one cannot at present identify 
a natural reason for e# to increase during the early afterglow phase, but future theoretical 
advances (see, e.g., Medvedev et al. 2005) might make it possible to study the evolution of 
shock-generated magnetic fields over time scales that are long enough to address this ques- 
tion. The value of e e might also be linked to the changing shock parameters (in particular, 
the shock Lorentz factor). However, in this case, again, our current level of understanding 
does not allow us to make a specific prediction. 

The possible increase of the afterglow radiation efficiency with time during the early 
phases of the X-ray light curve may also help to lower the estimate of the 7-ray radiative 
efficiency and thereby alleviate the constraints on the internal-shocks model. If at early 
times (t < t*) ex increased with t, then (by eq. [15]) -Ek,iso,o would be underestimated, 
and therefore the parameter / and the value of e 7 would be overestimated. The prompt 
emission from internal shocks could in principle be observable even if the very early radiation 
from the external (forward) shock is weak because of a low value of or of a possible 
suppression of e e when the shock Lorentz factor is high. This is because the value of e# in 
the internal shocks might be determined by a comparatively strong magnetic field advected 
from the central source (e.g., Spruit et al. 2001) and because (in contrast to the forward 
shock at the deceleration time, whose Lorentz factor is > 10 2 ) the internal shocks are only 
mildly relativistic. If the afterglow emission efficiency recovers its commonly assumed value 
(equation [4]) at t > t D rcak,2 then one could in principle have / ~ 1, k ~ 1, and e 7 ~ 0.5. 
In this case -Ek,iso remains constant while the flatfish flux decay is caused by an increase 
of ex with time. However, if even at late times the afterglow efficiency is only a fraction 
5 <C 1 of its standard value [where 5 might, for example, correspond to the fraction £ e of the 
electrons behind the forward shock that are accelerated to relativistic energies, which could 
conceivably be <C 1 (Eichler & Waxman 2005)] then k ~ 5/t s tandard ~ S 1 and e 7 ~ 5 (for 
^obs, grb ~ 1). In this case -Ek,iso is again constant with time and the flatfish decay phase is 
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due to ex increasing with time; however, the value of ex at late times is S times its standard 
value (given by equation [4]) and the true kinetic energy in the afterglow shock is a factor 
5 1 3> 1 larger than the usual estimate of ~ 10 51 erg. 

An alternative to the explanation of the flatfish segment of the light curve in terms of 
a prolonged energy injection into the forward shock of a basic jet model is the possibility 
(in a generalized jet model) of distinct spatial components, some of which only contribute 
to the afterglow emission at later times. This situation could arise if (1) some of the ejecta 
have lower initial Lorentz factors that result in longer deceleration times, so only a small 
fraction of the injected energy is transferred to the shocked external medium early on, or if 
(2) radiation from components that do not move along our line of sight is strongly beamed 
away from us at early times, becoming visible only later on when deceleration causes the 
respective beaming cones to widen. Case (1) is exemplified by the two-component jet model 
considered in the next section, whereas case (2) might be realized in the "patchy shell" model 
of GRB sources (Kumar & Piran 2000) and in the "anisotropic afterglow efficiency" scenario 
outlined by Eichler & Granot (2005). 9 



4. The Two-Component Jet Model: A Case Study 

As pointed out in § 1, the two-component jet model was originally invoked by Peng et 
al. (2005) as a possible way of alleviating the pie-Swift constraints on the 7-ray emission 
efficiency. In this section we interpret the Swift results in the context of this scenario, using 
again the parameters k and / introduced in § 2 and affirming some of the conclusions reached 
in that section. In this discussion we assume that the standard afterglow theory applies and 
that the microphysical parameters do not change with time. This formulation is used to 
demonstrate that the recent observations are inconsistent with parameter values of the two- 
component jet model that could lead to a lower inferred magnitude for e 7 , reinforcing the 
result that, in the context of the standard afterglow theory, the Swift measurements have 
tightened the constraints on the prompt-emission efficiency We go on to show, however, 
that the two-component model can nevertheless account for the X-ray afterglow light curve 
of GRB sources, including the flatfish early-time segment. 

The generic two-component jet model consists of a narrow and initially highly relativistic 
(conical) outflow from which the prompt emission originates, and a moderately relativistic 
flow that decelerates at a significantly later time and occupies a wider (coaxial) cone. The 



9 Case (2) might also be relevant to the two-component jet model if the line of sight to the observer lies 
within the solid angle subtended by the narrow component but outside the inner edge of the wide component. 
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narrow component is the source of the prompt emission (which is observed as a GRB if the 
observer's line of sight lies within, or very close to, its opening solid angle), whereas the wide 
component only makes a contribution to the afterglow emission (which becomes observable 
after it decelerates). Letting Ei, E^, 9i, and 77, stand for the total initial energy, initial 
kinetic energy, opening half-angle, and Lorentz factor of the two components (with % — n, w 
corresponding to narrow and wide, respectively), we have 

2 

-^7,iso £7-^11, iso ^7^ n 7)2 ' (1^) 

n 

The parameter k defined in equation (3) can be expressed in this context by the relation 

E°%Y = e obs,GRBe 7 £n (j^j ~ « (E k , n + ^k,w) , (17) 

where E k ^ n = (1 — e 7 )_E n and E k:W ks E w . Equation (17) incorporates the fact that the kinetic 
energy responsible for the late afterglow emission is determined empirically by assuming 
a jet of half-opening angle 9 W , or equivalently that the isotropic-equivalent kinetic energy 
inferred from late time afterglow corresponds approximately to the total kinetic energy over 
the fraction of the total solid angle occupied by the wide component. 

The energy that determines the early afterglow phase is 

2 

-Eearly,iso E k ^ n —^ -El<:,n,iso i (18) 
n 

whereas the late afterglow phase is determined by 

2 

£late,iso = (-Ek.n + £k,w) • (19) 

The evidence from the Swift observations that we do not see the early afterglow emission 
above the rapidly decaying tail of the prompt emission, and that even when it shows up it 
is rather weak, implies that the early isotropic kinetic energy -E ea ri y ,iso is a factor / ~ 10 
smaller than the late isotropic kinetic energy -Ei a te,iso> where we reintroduced the parameter 
/ defined in equation (3). This implies that 



Ek,n E k „ + E kw E kw 



(20) 



01 foi foi 

which in turn implies that 

£k,n < £k,w • (21) 

Hence we can delete the term involving E^ n from the right-hand side of equation (17) and 
obtain 

^obs.GRB E^ w ( 9 n \ E k 

,w,iso 



= K ^w =.r™. (22) 

(1 — e 7 J _C/k,n \fw/ ^k,n,iso 
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Expressing now E^ w in terms of E^ n using this last equation and substituting into equation 
(20) (which becomes simply / fa E ktWtiso / E k ^ iso when one neglects E kn on the r.h.s.), we 
rediscover the first relation in equation (3), which can be expressed in the form of equation 
(13) to yield e 7 ~ 0.9 for / ~ 10 and (/t/e obSiG RB) ~ 1- 

Equation (22) with e 7 e bs,GRB ~ 1 implies that K,E k:W /9^ fa E n /9l, and hence, given that 
£k, w ~ E w and taking k to be ~ 1, that the narrow and wide jet components initially have 
comparable isotropic-equivalent energies. The ratio of true energies of the two components 
is initially E w /E n fa k,~ 1 (6 w /6 11 ) 2 ~ (4 — 9)/t~ 1 for reasonable ratios of the opening half- 
angles. This ratio further increases by a factor (1 — e 7 )~ 1 fa nf (see eq. [13]) during the 
prompt-emission phase. 

As argued by Peng et al. (2005), this model could reduce the inferred value of e 7 if 
E kjW /E k>n > 1 and £ , k i w i i SO /-Ek,n,iso < 1- However, the Swift results, as expressed by equation 
(20), demonstrate that the latter ratio is equal to / ~ 10, and hence that this possibility 
is not realized in practice. Peng et al. (2005) also discussed the ability of this scenario 
to account for certain "bumps" in the late-afterglow light curve of several pre-Swift GRB 
sources. We now show that this model can similarly account for the early-time flattening 
of the X-ray light curve. The required model parameters are, however, distinctly different 
from those considered by Peng et al. (2005); in fact, the fits that we obtain reinforce the 
conclusion that -Ek,w,iso/i?k,n,iso (ss /) must be >> 1. 

The Lorentz factor of the narrow jet component does not directly enter into our modeling 
of the light curve, although the usual "compactness" arguments for the prompt emission (e.g., 
Lithwick & Sari 2001) imply that its value is r] n > 10 2 . Our interpretation of the flatfish 
segment of the light curve in the context of this scenario is that it largely corresponds to the 
emergence of the wide component around its deceleration time tdcc.w (-Ek^iso/TW) 1 ^ 3 ^ 8 ^ 3 , 
where n ext is the particle density of the ambient medium at the deceleration radius. Owing 
to the strong dependence of tdec on 77, we can constrain the value of i] w within a factor of 2 
or so. 

Figure 1 shows a tentative fit to the X-ray light curve of GRB 050315 with the two- 
component jet model, whereas Figure 2 demonstrates the effect of modifying the model 
parameters. The extended flat segment of the light curve together with the requirement 
that the contribution from the narrow component at t/(l + z) ~ 200 s does not overproduce 
the observed flux imply / fa E K<lso / E KnM > 30. This suggests that / ~ f max ~ 30 (see 
Table 2) and that both / min and / max might have been somewhat underestimated for this 
event (since the fit to the data should produce / > / m i n , suggesting that / min ~ 30 for 
GRB 050315, which is higher than the value of / m i n =11 derived by Nousek et al. 2005 and 
shown in Table 2). For our fiducial parameter values (e e = 0.1, e# = 0.01, and p fa 2.0 — 2.1; 
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see Table 2 and Fig. 1) we find (using eq. [5]) k 0.034 (see Table 2). For this value of k the 
product Kf is ~ 1, which implies (from eq. [13] with e b s ,GRB ~ 1) that e 7 ~ 0.5. The fit to 
the data shown in Figure 1 incorporates an SSC component whose contribution in this case 
turns out to be comparable to that of the synchrotron component at t* = 10 hr, resulting 
in a decrease by a factor of ~ 2 in the estimate of £k,iso(10 hr) (and in a corresponding 
increase in the estimate of k) in comparison with the values inferred by using equation (5) 
(which only takes account of the synchrotron contribution). The actual numerical difference 
between these two estimates is, in fact, slightly larger (a factor ~ 2.5), reflecting the fact 
that the fit employs a numerical scheme that is not identical to equation (5). All in all, the 
model fit shown in Figure 1 yields / w £ k ,w,iso/ #k,n,iso = 30, k w £°^ /£ kjWiiso w 0.086, and 
K f w S^/Ek.n.fao ~ 2.6, which implies e 7 ss 0.72 for e obS)G RB = 1. 

One realization of a two-component relativistic outflow of the type considered here is 
an initially neutron- rich, hydromagnetically accelerated jet (Vlahakis et al. 2003). In this 
picture the narrow component consists of the originally injected protons that are gradually 
accelerated to r) n ; the injected neutrons decouple from the original proton component when 
the jet Lorentz factor reaches i] w and eventually decay to form a distinct (wide) proton 
component. In the illustrative solution presented in Vlahakis et al. (2003), i] w ps 15. The 
dashed curve in Figure 2 demonstrates that, to be consistent with the data, this value of 
r/ w requires a very high external density, n ext ~ 10 3 ' 5 cm -3 (although even for this value 
the fit is not as good as for the parameters adopted in Fig. 1). This can be understood 
from the parameter dependence of tdcc,w/(l + z), for which the model fit implies a value of 
~ 2 x 10 3 s. To reproduce the observed flux, the value of -Ek,w,iso cannot be too low. In 
fact, in the model fit represented by the dashed curve we have adopted the lowest possible 
value of this quantity, corresponding to an equipartition between the electron and magnetic 
field energy densities. With the values of £dcc,w and -Ek,w,iso thus fixed, the inferred external 
density becomes very sensitive to the value of the Lorentz factor (n ext oc i] w 8 ). Intermediate 
options, with a somewhat higher Lorentz factor (i] w = 21) are also shown in Figure 2, both 
for the uniform-density case used in the previous fits (dash-dotted curve) and for an r~ L5 
density profile (where r is the distance to the source; dotted curve). We interpret the break 
in the light curve of GRB 050315 at t/(l + z) ~ 1 day as the jet break time of the wide 
component, tj et;W /(l + z). Since the jet break time is particularly sensitive to the value of 
the opening half-angle, £j e t,w oc (-Ek^iso/'ftcxt) 1 ^ 3 ^ 3 , this allows us to constrain the value of 
9 W . Given that tj e t,w/£dec,w oc (f] w tW ) 8 ^ 3 , it is seen that any reduction in the fitted value of 
i] w requires an increase in 9 W by a similar factor. 

As was already noted in Peng et al. (2005), another realization of this type of an outflow 
is potentially provided by a relativistic, baryon-poor jet, which is driven electromagnetically 
along disk-anchored magnetic field lines that thread the horizon of a rotating black hole, 
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and which is "contaminated" by neutron diffusion from a baryon-rich disk wind. In the 
original version of this scenario, which was proposed by Levinson & Eichler (1993, 2003; see 
also van Putten & Levinson 2003) and recently studied numerically by McKinney (2005a,b), 
the narrow and wide components correspond to the baryon-poor and baryon-rich outflows, 
respectively. However, the predicted Lorentz factor of the disk wind is too low to be consistent 
with the rather high value of r] w inferred from our model fits. An alternative possibility is 
that the wide component corresponds to a hydromagnetically accelerated baryon-rich disk 
outflow of the type modeled by Vlahakis & Konigl (2003a,b), whereas the narrow component 
corresponds (as in the Levinson & Eichler picture) to a baryon-poor Blandford & Znajek 
(1977)-type outflow. 10 It should be noted that in either one of the above two-component jet 
realizations the 7-ray emitting component is associated with an initially Poynting-dominated 
outflow. This could in principle make it possible to account for the relatively high emission 
efficiency that the Swift results seem to imply (see § 5). 



5. Summary and Conclusions 

We have shown that the 7-ray efficiency implied by the Swift observations is model- 
dependent and can vary over a wide range (from typical values of e 7 ~ 0.9 or higher to 
e 7 ~ 0.1 or lower) depending on the adopted model assumptions. The 7-ray efficiency has 
been expressed in terms of observable quantities (see eqs. [3] and [13]), namely k and /, 
where k relates the 7-ray emission to the late-time afterglow emission (and was therefore 
available in the pie-Swift era) and / relates the early- and late-time afterglow energies (and 
therefore became available only with the launching Swift). We have shown that there is 
a large uncertainty in the values of both k and /, which translates into a corresponding 
uncertainty in the value of the 7-ray efficiency, e 7 . In the following discussion we make the 
conservative assumption that we observe most of the emitted energy in 7-rays (e b s ,GRB ~ 1); 
if a significant fraction of the total radiated energy is emitted outside of the observed photon 
energy range (e.g., at higher energies) then this would increase the required value of e 7 (see 
equation [13]). 

The kinetic energy of the afterglow shock has been estimated in the pre-Swift era using 
broad-band afterglow fits for a small number of GRB sources that had the best available 
broad-band afterglow data, yielding a typical value of k ~ 1 with a large scatter of almost 



10 It is in principle also conceivable that the baryon-poor outflow could develop an internal structure that 
would correspond to the two outflow components considered here, but there are at present no quantitative 
results to support this conjecture. 
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an order of magnitude (see § 2.3). Substituting the values of the microphysical parameters 
(e e , €b and p) that were derived from these fits into our equations for E ktiso generally yielded 
somewhat lower values of k. Finally, using our equations with fiducial values of the micro- 
physical parameters (e e = 0.1, e# = 0.01, and p = 2.2) gives a typical value of k ~ 0.1 — 0.2, 
both for the pie-Swift and the Swift GRB samples that we use, with a moderate scatter. 
Specifically, for our pie-Swift (Swift) sample, (log 10 k) = —0.75 (—0.82) corresponding to 
k = 0.18 (0.15) and o"i ogl() K = 0.60 (0.63). Obviously, the choice of fiducial values for the 
microphysical parameters is somewhat arbitrary and affects the resulting value of k. Higher 
values of the microphysical parameters e e and e# (e.g., e e ~ 0.3, e# m 0.08) are required 
in order to obtain an average value of log/t, using our equations, similar to that derived 
from pie-Swift broad-band afterglow fits. Altogether, there is almost an order-of-magnitude 
uncertainty in the typical value of k (which ranges from ~ 0.1 — 0.2 to ~ 1). 

Even if we adopt the high typical values of k (~ 1) inferred from pre- Swift afterglow 
fits, it is important to keep in mind that these values have been estimated on the basis of 
the standard assumptions of afterglow theory. Changing these assumptions would modify 
the inferred value of k for the same fits. For example, as pointed out by Eichler & Waxman 
(2005), if only a fraction £ e < 1 of the electrons are accelerated to relativistic energies, 
then there is a degeneracy where the same observable quantities are obtained for e e — > £ e e e , 
e_e — ► ie^B, n cxt — > CT^ext, and -Ek,iso — > CT^k.iso- Since this increases the inferred value of 
^k,iso by a factor of k is reduced by the same factor in comparison with the estimate 
from the standard theory (which uses £ e = 1). 

An alternative way of reducing the inferred value of e 7 was proposed by Peng et al. 
(2005) in the context of the pie-Swift observations. Specifically, they considered a two- 
componet outflow model with parameters that effectively corresponded to the parameter / 
having a value < 1. This parameter choice appears to be inconsistent with Swiff s detection 
of an early flatfish decay phase in the X-ray light curve, which, when interpreted in the 
context of the standard afterglow theory as arising from a gradual increase with time of 
E^iso, typically implies / ~ 10 (and, in some cases, values of / that are as high as ~ 
10 2 — 10 3 ). It is worth noting in this connection that the existence of a two-component GRB 
jet model can be plausibly expected on various theoretical grounds and has been suggested 
independently on the basis of fits to several pie-Swift afterglows (see discussion in Peng et 
al. 2005). Furthermore, the Swift observations by no means rule out this model, although 
they can be used to constrain its parameters. We have illustrated this fact through the fit 
to the X-ray light curve of GRB 050315 that we performed in § 4 within the framework 
of this model. This fit has yielded plausible ranges for the kinetic energies and opening 
angles of the two components as well for as the Lorentz factor of the dominant (wide) 
component. A key conclusion from this fit is that the kinetic energy of the wide component 
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is much larger than that of the narrow one (£k,w/-Ek,n ~ 10 2 ). Physically, the narrow and 
wide components could conceivably correspond to a baryon-poor, black-hole-driven outflow 
(Levinson & Eichler 1993, 2003) and a baryon-rich, disk-driven outflow (Vlahakis & Konigl 
2003a,b), respectively, although this remains to be demonstrated. We also note that the two- 
component jet parameters derived in § 4 were based on standard assumptions; they could 
change if the underlying assumptions (involving, for example, the values and time constancy 
of the microphysical parameters) were altered. 

In conjunction with the k ~ 1 estimates of the standard broad-band afterglow fits, 
values of / > 10 imply 7-ray radiative efficiencies e 7 > 0.9. Such high efficiencies would 
be essentially impossible to achieve in any scheme, such as the internal- shocks model, that 
tapped the kinetic energy of the outflow for radiation. An alternative possibility that has 
been discussed in the literature is the direct transfer of Poynting flux (which evidently is also 
a major contributor to the flow acceleration — e.g., Drenkhahn & Spruit 2002; Vlahakis & 
Konigl 2003a) into nonthermal radiation (e.g., Usov 1994; Thompson 1994). It is at present 
unclear how to assess the efficiency of this process. There are two generic possibilities: 
dissipative fronts and magnetic reconnection sites. The first option corresponds to overtaking 
collisions of magnetically dominated relativistic streams and is not expected to result in high 
radiative efficiencies (e.g., Romanova & Lovelace 1997; Levinson & van Putten 1997). The 
second case would require magnetic field orientation reversals and would most naturally arise 
in a pulsar-type outflow from a rapidly rotating neutron star (e.g., Spruit et al. 2001). In this 
case it is, however, still unclear whether radiative efficiencies > 0.9 could be attained even 
under the most favorable assumptions about the field reconnection rate (Drenkhahn & Spruit 
2002), and it has in fact been suggested that the reconnection rate might be self-limiting 
(Lyubarsky & Kirk 2001). 

An "intermediate" situation could prevail if / ~ 10 but k ~ 0.1, reflecting the pos- 
sibility that k was overestimated by the pie-Swift afterglow fits (perhaps because some of 
the assumptions of the standard theory do not hold — e.g., £ e ~ 0.1 rather than £ e = 1). 
Alternatively, k could be ~ 1 but f — 1, corresponding to the early flatfish decay phase 
reflecting an increase with time of the X-ray afterglow efficiency ex (due, e.g., to p being < 2 
or to an increase with time of e e or e#; see § 3) rather than an early increase in E^ iso . In 
either one of these cases the inferred 7-ray radiative efficiency would be reduced to e 7 ~ 0.5. 
Although this value is less extreme than the estimate discussed in the preceding paragraph, 
it is worth noting that it is still fairly restrictive for the internal-shocks model, in which it 
could potentially be attained only if all of the following conditions (already summarized in 
Peng et al. 2005) are satisfied: (1) the ratio between the maximum and minimum initial 
Lorentz factors of the ejected shells is large enough (> 10); (2) the distribution of initial 
Lorentz factors is sufficiently nonuniform; (3) the shells are approximately of equal mass 
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and their number is large enough (> 30), and (4) the fraction of the dissipated energy that 
is deposited in electrons and then radiated away is sufficiently high (e ei GRB e rad,GRB 
with a similar constraint applying to the fraction e bs,GRB of the radiated energy that is 
emitted as the observed 7-rays (see Beloborodov 2000 and Kobayashi & Sari 2001). Only if 
both f ~ 1 and k ~ 0.1 were satisfied (which could occur, for example, if e e or e# initially 
increased with time and £ e were ~ 0.1) would the inferred value of e 7 drop to ~ 0.1 and be 
compatible with the values that are expected to arise under less constrained circumstances 
in the internal-shocks model. 
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Table 1. Estimates of k for pre- Swift GRBs with known redshifts 
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Note. — Estimates for i?k,iso(10 hr) = i?k,iso,i0hr and k = i?° b ;| o /-E ; k,iso,l0hr for the GRBs with known redshift from the Bloom et al. (2003) and Berger et al. (2003) 
samples, t L x>isOjl0hr = L x , iso (10 hr) from Table 2 of Berger et al. (2003), * B* s so in the 20 - 2000 kcV range form Table 2 of Bloom et al. (2003), H calculated using eq. 
(5) with our fiducial values of the microphysical parameters: e e = 0.1, eg = 0.01, and p = 2.2. * substituting the best fit values of the microphysical parameters from 
Table 2 of Panaitescu & Kumar (2002) into eq. (5), * the value for -Ek,iao(10 hr) is taken to be fa 0.5-E]j,iso,o> where -Ek.iso.o i s taken from Panaitescu & Kumar (2002) 
and the factor of 0.5 accounts for the average factor by which the energy is reduced relative to -Ek.iso.o due to radiative losses (A. Panaitescu, personal communication), 
* these values are for a fit to a uniform external density; PK02 get a significantly better fit to a wind density for which we derive n = 0.47 and 0.17 instead of 0.42 and 
0.14, respectively, * in these cases Panaitescu & Kumar (2002) find p < 2 which introduces an uncertainty through the extrapolation that is involved in the expression 
for the numerical coefficient in eq. (9); for GRB 991216 PK02 find that the X-ray emission is from electrons with 7 e > 7 max where dN/dje <x 7 e q with p < 2 < q so 
that we cannot readily substitute their results into our equations, ^ substituting the best fit values of the microphysical parameters from Table 1 of Yost et al. (2003) 
into eq. (5), <> £ k ,iso,ld = #k,iso(l day) = E° b ?J Kld from Table 1 of Yost et al. (2003). <><> K ld = E° b ?JE kiiBoM . 
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Table 2. Estimates of / and k for Swift GRBs with known redshifts 
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Note. - - The estimates for / = E KisoAOhr /E^ sofi and k = E° h *J E KisoAOhr for the GRBs 
with known redshift from the Nousek et al. (2005) sample; t Lx,iso,iohr = -^x,iso(10 hr) in the 
2 — 10 keV range at (both the time and the photon energies measured in the cosmological 
frame of the GRB) from Table 2 of Nousek et al. (2005), t £°^ Q in the 20 - 2000 keV range 
(in the GRB's cosmological frame) form Table 2 of Nousek et al. (2005), ™ calculated using 
eq. (5) with e e = 0.1, e# = 0.01, and the values of p from this Table, * estimated using the 
measured spectral slope in the X-rays (Nousek et al. 2005) and attempting to fit it into the 
range 2 < p < 3 if allowed within the errors on the spectral slope, * k = i?°^ /£'k i iso,iohr 
estimated using the values from this Table, 9 / min and / max are taken from Table 3 of 
Nousek et al. (2005), and are estimated using the measured X-ray flux at tbrcak.i) an d the 
extrapolated X-ray flux at T G rb, respectively (see Nousek et al. 2005 for details). 
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t / (1 +z) in seconds 

Fig. 1— Tentative fit to the X-ray light curve of GRB 050315 (from Nousek et al. 2005) 
with the two- component jet model. The numerical code used to calculate the light curve is 
essentially Model 1 of Granot & Kumar (2003), which neglects the lateral spreading of the 
jet and includes synchrotron self-Compton (SSC) emission. In addition to the total light 
curve (thick solid line) we also show the separate contributions of the different components: 
the tail of the prompt emission (oc t~ 5 ), the narrow outflow, and the wide outflow. Here 
Ek = ^k,w + Ek,n is the total kinetic energy of the two components. The narrow and wide 
components occupy the non-overlapping ranges 9 < 9 n and 9 n < 9 < 9 W , respectively, in the 
polar angle 9 (measured from the symmetry axis); 9 ^ s is the viewing angle with respect to 
this axis. 
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Fig. 2. — The effects of varying the two-component jet model parameters with respect to 
those of the fit shown in Fig. 1 (reproduced in this figure by the solid line). Here E 5 i = 
E k /(10 51 erg) (with E k = E ktW + E k>n ), whereas E kjW / n denotes the ratio E k)W /E k>n of the 
kinetic energies of the wide and narrow components. 



